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In this talk, after reviewing the dynamical gluon mass generation mechanism within the pinch 
technique framework and its phenomenological predictions, we will introduce the modern formu- 
lation of the pinch technique which makes extensive use of the Batalin-Vilkovisky quantization 
formalism. In this framework a certain auxiliary function A fiv (q)- and its associated form factors 
G{q 2 ) and L(q 2 ) - play a prominent role. After showing that in the (background) Landau gauge 
Ajiv(<7) fully constrains the QCD ghost sector, we show that G(q 2 ) coincides with the Kugo- 
Ojima function u{q 2 ), whose infrared behavior has traditionally served as the standard criterion 
for the realization of the Kugo-Ojima confinement mechanism. The determination of the behavior 
of G(q 2 ) (and therefore of the Kugo-Ojima function) for all momenta through a combination of 
the available lattice data on the gluon and ghost propagators as well as the dynamical equation 
G(q 2 ) satisfies, will be then discussed. In particular we will show that (z) in the deep infrared the 
function deviates considerably from the value associated with the realization of the Kugo-Ojima 
confinement scenario, and (//) establish the dependence on the renormalization point of u(q 2 ), 
and especially of its value at q 2 = 0. 
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In the last two years, ab-initio lattice gauge theory computations using extremely large vol- 
umes have firmly established that (in the Landau gauge) the QCD gluon propagator and the ghost 
dressing function are infrared (IR) finite and non- vanishing Specifically choosing an/?^ type 

of gauge and defining the gluon propagator cofactor A, and the ghost dressing function F as 

a mv(<?) = -i 

where Pj lv (q) =g^ v — q^qv/q 2 is the transverse projector, A~ l (q 2 ) = q 2 + ill(q 2 ) [with Yl llv {q 2 ) = 
P tlv {q)Yl{q ) the gluon self-energy], and D{q 2 ) is the ghost propagator, lattice results tells us that 
(Euclidean space) 

A _1 (0)>0, and F(0)>0. (2) 

The issue of explaining these clean lattice results from the point of view of the continuum formu- 
lation of the theory has therefore become an increasingly interesting topic, for obtaining them is 
bound to expose a QCD fundamental dynamical mechanism at work. 

Indeed, such a mechanism is provided by the dynamical generation of a gluon mass [||] 
through the non-perturbative realization of the well-known mechanism described long ago by 
Schwinger Schwinger mechanism shows that if for some reason the dimensionless vacuum 
polarization Yl(q 2 )/q 2 behaves as a simple pole with positive residue /I 2 at q 2 = (an there is 
no physical principle precluding this possibility) then A~'(g 2 ) = q 2 + /I 2 : Thus the vector meson 
(which is massless in the absence of interactions) becomes massive, with A _1 (0) = ft 2 . When 
the theory is strongly coupled, as it happens with QCD in the IR, strong binding may generate 
zero-mass bound state excitations which, notwithstanding the fact that they do not generate from 
the spontaneous breakdown of any symmetry, acts like massless, composite and longitudinally 
coupled (dynamical) Nambu-Goldstone bosons [j|]. 

The implementation of this mechanism within the pinch technique (PT) framework |^ [7||, 
gives rise to two complementary effects, which appears at the level of the QCD Schwinger-Dyson 
equations (SDEs), and of the effective low-energy theory which describes QCD in the IR sector, 
eventually providing a confinement mechanism. 

• Schwinger-Dyson equations. The systematic exploitation of the underlying BRST symmetry 
provided by the PT originate drastic modifications to the Green's functions of the theory 
and the corresponding SDEs which describes their dynamics; in particular the new SDE 
obtained for the gluon propagator lend itself to a novel truncation scheme that respects gauge 
invariance at every level of the dressed loop expansion [Q]. Once this scheme is used together 
with the assumption that the three-gluon vertex contains dynamical massless poles (~ 1 /q ) 
triggering the Schwinger-mechanism, one obtains IR finite solutions for the propagator 
These solutions are of the type A -1 (q 2 ) = q 2 + m 2 {q 2 ) with the mass depending non-trivially 
on the momentum transfer, and with (obviously) A _1 (0) = m 2 (0) = m 2 , > (Fig. [T], left 
panel). In addition of taming the IR divergences intrinsic to perturbation theory (Landau 
pole) this mass forces F(q 2 ) to stay IR finite with F(0) > (Fig. |], right panel), and no 
enhancement is found [|l0|]. 

• Effective low-energy theory. An effective low-energy field theory for describing the gluon 
mass is the gauged non-linear sigma model known as "massive gauge-invariant Yang- 
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Figure 1: Left Panel: The numerical solution for the gluon propagator from the PT modified SDE (black 
solid line) compared to the lattice data of [01, Right panel: The ghost dressing function p 2 D(p 2 ) obtained 
from the SDE. In the inset we show the lattice data for the same quantity; notice the absence of any enhance- 
ment in both cases. 



Mills" Qll|], with Lagrangian density 

1 



^M = ^ v -mlTv[A ll -g- l U{B)d pL U-\e)\ 



(3) 



where = ^ Y,a ^a^u> the X a are the SU(3) generators (with TrX a Xb = 28 c ,b), and the N xN 
unitary matrix U (6) = exp [i^A a f3°] describes the scalar fields 6 a . Note that Jz? MYM is locally 
gauge-invariant under the combined gauge transformation 



= VA^V 



-l 



[dpV] V 



-l 



u' = u(e')=vu(6) 



(4) 



for any group matrix V = exp [/^A a ft)"(^)] , where a> a (x) are the group parameters. One 
might think that, by employing @, the fields B a can always be transformed to zero, but 
this is not so if the 6 a contain vortices. To use the Jz? M ym in ©, one solves the equations of 
motion for U in terms of the gauge potentials and substitutes the result in the equations for the 
gauge potential. One then finds the Goldstone-like massless modes mentioned above. This 
model admits vortex solutions [pi]], with a long-range pure gauge term in their potentials, 
which endows them with a topological quantum number corresponding to the center of the 
gauge group [Z N for SU(N)], and is, in turn, responsible for quark confinement and gluon 



screening. Specifically, center vortices of thickness ~ Mq , form a condensate because their 
entropy (per unit size) is larger than their action. This condensation furnishes an area law 
to the fundamental representation Wilson loop, thus confining quarks [0, 11]. In addition, 



the adjoint potential shows a roughly linear regime followed by string breaking when the 



potential energy is about 2mo, corresponding to gluon screening [12]. 



Thus, summarizing, in this picture the non-perturbative QCD dynamics generate an effective, 
momentum dependent mass, without affecting the local SU (3) invariance, which remains intact. 
This provides in turn very definite predictions about the IR behavior of the theory (all in agreement 
with lattice studies): (i) the gluon propagator is IR finite; (ii) in the Landau gauge the ghost remains 
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massless, but with a finite dressing function; the presence of center vortex solutions providing (Hi) 
an area law for the Wilson loop and (iv) a roughly linear behavior for the adjoint potential followed 
by string breaking. 

A different set of predictions is obtained within the Kugo-Ojima (KO) scenario, which also 
establishes a highly non-trivial link between confinement and the infrared behavior of some fun- 
damental Green's functions of QCD |[T3|]. In the KO confinement picture one starts by observing 
that the equation of motion for the gauge field can be written in the Maxwell-like form d y F£ v + 
{Gbrst, (@nc) a } = gJ 1 ^, with /" the Noether current of the global color symmetry, ^ the usual co- 
variant derivative, and <2brst the BRST charge operator. On the other hand, if one adds to a Noether 
current the derivative of an antisymmetric rank 2 tensor (say e^/j^y]) tne resulting current is still 
preserved and the corresponding charge correctly generates rotations in color space. Thus one is, 
at least in principle, allowed to define the color charge Q" in the BRST exact form 

Q a = J d 3 x (j$ - l ~d V Fov) = J d 3 ^{e BR sT, {%c) a }. (5) 

Confinement is then a direct consequence of the above relation, since for any physical state speci- 
fied by the condition <2 B RST|phys) = one has (phys|<2 a |phys'} = 0, which implies that all physical 
states are color-singlets fll3|]. 

The problem is however that the volume integral in Eq. ^ does not converge, due to the 
presence of massless one particle contributions to d v F^ v , and {<2 B rst, (^u c ) a } (the so-called 
quartet mechanism). Without entering into any detail, one has that a solution to this problem 
is provided by introducing some suitable weights v, w and u respectively, so that a well defined 
charge is given by 

Q a = J d 3 x (jg + ^d v Fg v ) , gv = -w + (1 + u) (6) 

where the relation represents the condition between the different weights for the cancellation of the 
aforementioned one-particle massless contributions. Thus, requiring that Eq. (H) coincides with the 



BRST exact expression (ph implies v/w = — l/g, and therefore the KO confinement criterion [13] 



1 + u = 0. (7) 

It turns out that in the Landau gauge the so-called KO parameter u is linked to the IR limit 
q 2 — > of a certain Green's function; more precisely one has 

j d 4 *e-^->'> (T [ {%c) n x (%c) n y j) = -^S^+P^q^utf), (8) 
lim u(q 2 ) = u(0) = u. (9) 

In addition, there is a powerful BRST identity relating the KO-f unction u(q 2 ) and the ghost dressing 
function F(q 2 ), namely [14, |I3| ] 

F- 1 (^ 2 ) = l + ^ 2 )+w(^ 2 ), (10) 
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where for now w(q 2 ) is an unspecified function, which under very general conditions is such that 
w(0) = 0. Then, Eq. (0) tells us that the KO confinement scenario predicts an IR divergent ghost 
dressing function, at odds with the aforementioned large volumes lattice simulations 1 . 

An issue that to the best of our knowledge has never been thoroughly addressed in the KO sce- 
nario is how renormalization affects the proof of the central identity (^). On the other hand, and at 
a less formal level, one should notice that once the relations (|9|) are proved, the KO function is on a 
par with any other QCD Green's function: quantum corrections will set in, and the whole procedure 
of regularization and renormalization should be applied. This in general implies the unavoidable 
appearance of a dependence in the Green's function (^) on a "sliding scale" ju. This dependence 
disappears only (z) when combining individual Green's functions to form observables, such as S 
matrix elements, or (z'z) when forming very special (and very well studied) products of Green's 
functions and the gauge coupling of the theory. The latter are the so-called renormalization-group 
(RG) invariant combinations, with the product e A in QED constituting probably the most cele- 
brated text-book case. Whether a product of Green's functions forms a RG-invariant combination 
or not is determined on formal grounds, from the Ward-Takahashi or Slavnov-Taylor identities 
satisfied by the quantities involved {e.g., the famous Z\ = Z2 of QED). The Green's function of 
Eq. @ is definitely not an RG-invariant []l7|], and therefore picks up a non-trivial dependence on 
fi. Thus, one has u(q 2 ,fL 2 ), and in particular, in the deep IR limit, u = u(0,fx 2 ). 

Now, if the KO confinement criterion (|7]) were satisfied, the ghost dressing function must 
diverge as q 2 — > 0, due to the BRST identity (|10|). This will in turn make the ji. dependence of the 
KO function (which is however still there) irrelevant, since u(q 2 ,fL 2 ) will be then driven to -1 in 
the deep IR limit for any value of fL. On the other hand, lattice simulations tells us that the ghost 
dressing function is finite, and it is interesting to determine explicitly the /I -dependence of the KO 
function. 

Quite remarkably this issue can be thoroughly studied within the modern formulation of the 
PT (by means of the Batalin-Vilkoviski quantization formalism [18]) exploiting in particular the 
PT correspondence [Q, [HJ] with the background field method (BFM) [20]. Indeed, it turns out 
that [|l^, [l7|], in the (background) Landau gauge, the Kugo-Ojima function coincides with the 
form factor G(q 2 ) multiplying g^ v in the Lorentz decomposition of a certain auxiliary function 
A[i V (q) which enters in all the so-called "background-quantum" identities [21], i.e., the infinite 
tower of non-trivial relations connecting the BFM Green's functions to the conventional ones (e.g., 
calculated in the gauges). In addition, G(q 2 ) is a key element in the aforementioned new SDEs 
that can be truncated in a manifestly gauge invariant way. 

Let us recall that in the BV formulation of Yang-Mills theories [|l^], one starts by introducing 
certain sources (called anti-fields in what follows) that describe the renormalization of compos- 
ite operators; the latter class of operator is in fact bound to appear in such theories due to the 
non-linearity of the BRST transformation of the elementary fields. In much the same way, the 



1 The same prediction is obtained when implementing the (original) Gribov-Zwanziger (GZ) horizon condition [ |li| ] : 
in the IR region the ghost propagator diverges more rapidly than at tree-level. Furthermore, it has been also argued that 
the Landau gauge gluon propagator should vanish in the same limit. The prediction for the gluon propagator in the 
KO scenario is instead that it is less divergent than the tree-level expression which, evidently, encompasses the IR-finite 
gluon propagator as a special case, even though, up until recently, the focus had been placed rather on the "vanishing" 
solutions, given that they satisfy simultaneously both the KO and GZ requirements. 
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Figure 2: Diagrammatic representation of the functions A and H. 



quantization of the theory in a background field type of gauge requires, in addition to the afore- 
mentioned anti-fields, the introduction of new sources which couple to the BRST variation of the 



background fields pi]. These sources are sufficient for implementing the full set of symmetries of 
a non-Abelian theory at the quantum level, and in the case of quarkless SU (N) QCD, lead to the 
master equation 



/ 



d 4 x 



sr sr dr dr 



<5A* m SA™ 



8c m 8c 



— + a m 
8c™ + " 



( 8T 


8r y 




8AJl)_ 



(11) 



In the formula above, T is the effective action, A* and c* the gluon and ghost anti-fields, A is 
the gluon background field, and Q. the corresponding background source; finally B denotes the 
Nakanishi-Lautrup multiplier for the gauge fixing condition. 

To determine the complete algebraic structure of the theory we need two additional equations. 
The first one is the Faddeev-Popov equation, that controls the result of the contraction of an anti- 
field leg with the corresponding momenta. In position space, it reads 



sr 

8c m 



sr 

8AI 



(12) 



where (W$>) m = d»® m + gf mnr A n ^ r [in the case of {W$>) m replace the gluon field A with a 
background gluon field A]. The second one is the anti-ghost equation formulated in the background 



field Landau gauge, which reads Jl5[ ] 
8T 



8c" 



8T 
8£lu 



(^A* ti ) m -f mnr c* n c r + f 



. ST . 
8B» l 



0, 



(13) 



This equation fully constrains the dynamics of the ghost field c, and implies that the latter will not 
get an independent renormalization constant. The local form of the anti-ghost equation ( |I3| ) is only 
valid when choosing the background Landau gauge condition (Si^Ap)" 1 = 0; in the usual Landau 
gauge, d^A'jl = 0, an integrated version of this equation is available. In fact, even though the results 
that follow will be derived for convenience in the background Landau gauge, they are valid also in 
the conventional Landau gauge of the . 
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Now, differentiation of the functional (11) with respect to a combination of fields containing at 
least one ghost field or two ghost fields and one anti-field (and setting the fields and sources to zero 
afterwards) will provide the Slavnov-Taylor identities of the theory. Differentiation with respect to 
a background source and background or quantum fields will provide, instead, the aforementioned 
background-quantum identities. Finally, differentiation of (12) and ([O]) with respect to fields and 



anti-fields or background sources give rise to relation among the different auxiliary ghost functions 
appearing in the theory. 

The important point is that, when carrying out these differentiations, the following function 
appears (Fig. ||) 

iA MV (<7) = r flflA .( 9 ) = g 2 C A [ H$D(k + q)AP°(k)H av (k,q), 

= i g ^G(q 2 )+i^L(q 2 ), (14) 

In the equations above, the color factor 8" 1 " has been factored out (as always in what follows), 
Ca represents the Casimir eigenvalue of the adjoint representation [Ca = N for SU(N)], and 
f k = jJ. 2e (2n)~ d Jd d k, with d = 4 — e the dimension of space-time. Finally, the function H^ v (k,q) 
(see Fig. ^ again) is in fact a familiar object, for it appears in the all-order Slavnov-Taylor identity 
satisfied by the standard three-gluon vertex; it is also related to the full gluon-ghost vertex T^(k,q) 
by the identity q v H^ v (k,q) = —iT^(k,q) [at tree-level, Hjfy = ig^ v and T^\k,q) = —q^,]. Indeed 
one finds the following results 



(i) When differentiating the functional ( |1 1[ ) with respect to a background source and a back- 
ground gluon, on the one hand, and a background source and a quantum gluon, on the other, 
we can combine the resulting equation and trade the resulting two point functions for the 
corresponding propagators to get the important background-quantum identity 

A- l =[l+G(q 2 )] 2 A- l (q 2 ). (15) 

The quantity A(q 2 ) appearing on the left-hand side of the above equation captures the running 
of the QCD p function, exactly as it happens with the QED vacuum polarization; for every 
value of the (quantum) gauge-fixing parameter one has (at the one loop level) A~ l {q ) = 
q 2 [l +bg 2 \og(q 2 / fx 2 )] where b = 11Ca/487T 2 . It is the identity (§J) that plays a central role 
in the derivation of the new set of SDEs that can be truncated in manifestly gauge invariant 
way. 

(ii) If we consider the background Landau gauge, differentiating the ghost equation ( fl2| ) with 
respect to a ghost field and a background source, and the anti-ghost equation ( |I3| ) with respect 
to a gluon anti-field and an anti-ghost, we get the relations 

Tcc{q) = -iq V r cA * (q) Tga^ (q) =q^+ q v A^ v (q), (16) 

r c A* v (?) = qv + q^A^viq), F cd (q) = -iq^T^ (q). (17) 



Next, contracting the first equation in (17) with q v , and making use of the first equation 



in (|T6|), we see that the dynamics of the ghost sector is entirely captured by A^ v (q), for one 
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Figure 3: Connected components contributing to the function G"J v (q). 



has iT C c(q) = q 2 + q^q v A^ v (q). Then, introducing the Lorentz decompositions F C A^(q) 
q jl C(q 2 ) and (q) = q^E(q 2 ) we find the identities [|T|] 



C(q 2 )=E(q 2 )=F-\q 2 



F-\q 2 ) = \+G{q 2 )+L{q 2 ). 



(18) 



Recalling that the dimension of the gluon anti-field A* is three, while the dimension of the 



Q, source is one, power counting shows that all functions appearing in Eqs. (16) and (17) 
are divergent, and in particular that the divergent part of A llv (q) can be proportional to g^ v 
only [jl5L p2h, so that L{q 2 ) is ultraviolet finite. 

(in) In the background Landau gauge the function appearing on the lhs of Eq. ([|) is precisely 
given by 

8 2 W 

where W is the generator of the connected Green's functions, and the two connected diagrams 
contributing to are shown in Fig. |3[ Factoring out the color structure and making use of 
the identities (|T8b one has 



-i% v (q) = A^(q)+r a ,c(q)D(q 2 )r A * c (q) 



" q 2 



+ P llv {q)G{q 2 



(20) 



Passing to the Euclidean formulation, and comparing with Eq. <jsj> , we then arrive at the 
important equality 



u(q ) = G(q 



(21) 



Then, the usual KO confinement criterion may be equivalently cast in the form: 1 + G(0) = 0; 
moreover we see that the unspecified function w(q 2 ) appearing in Eq. (|l0|) coincides in fact 
with the L(q 2 ) form factor appearing in (|l4]). 

We thus see that in the (background) Landau gauge the single function G(q 2 ) is an extremely 
interesting object to study, since its IR behavior determines much of the IR behavior of the theory. 
Now, keeping the vertices T ll (k,q) and H^y^k^q) at their tree-level values, one finds that the dy- 
namical equations for F(q 2 ), G(q 2 ) = u(q 2 ) and L(q 2 ) are completely determined by the gluon and 



ghost propagator [22] 



F-\q 2 ) = l+g 2 C A j k 



(k-q) 



21 



L(q 2 



g 2 c A 



g 2 c A 



k 

2 + 



1-4 



k 2 q 2 

(k-q) 2 ' 
k 2 q 2 

(k-q) 



k 2 q 2 



A(k)D(k + q), 
A(k)D(k + q), 
A(k)D(k + q). 



(22) 
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u = 3.0 GeV and a(n ) = 0.388 
u = 3.6 GeV and a(n 2 ) = 0.330 
f. = 4.3 GeV and o(|i Z ) = 0.295 




q'[GeV] 




q'[Ge\T] 



Figure 4: Left panel: —u(q 2 ) determined from Eq. (24) at different renormalization points fi = 3.0, 3.6,4.3 
GeV. Right panel: Same as in the previous panel but this time for L(q 2 ). 



Now, as discussed in detail in [p2|], these (unrenormalized) equations must be properly renor- 



malized, i.e., in such a way as to preserve the validity of the BRST identity in (|10[), which should 
not be deformed by the renormalization process. Note in fact that Eq. ( JTo| ) constrains the cutoff- 
dependence of the unrenormalized quantities involved 2 . Specifically, if we denote by Z c the ghost 
wave-function renormalization constant, with Z c Fq 1 = F~ l and with Z\ the (yet unspecified) 
renormalization constant of the function A flv (q), with 3 Zp\g^ v + Aq V ] = g^ v + A^ v , one finds 
that in order to preserve the identity (10) one has to impose Za = Z c [g2|]; as a result, one finds the 
relation 

Z C (A 2 „U 2 )[1 +u {q 2 ,A 2 ) +L Q {q 2 ,A 2 )} = \+u{q 2 ^ 2 ) +L{q 2 ^ 2 ). (23) 

Imposing then the renormalization condition F(/j. 2 ) = 1, going to Euclidean space, setting q 2 = 
x, k 2 = y and a s = g 2 /4n, and implementing the standard angular approximation, one finds the 
renormalized equations 



F-\x) 
1 + u(x) 
L(x) 



167T 
167T 



F{X) J% y (3 - Z) A(y) + J"dy U - ^ A(y)F(y) 



x 

Fix 



Yin 



Fix) 



rx poo 

/ dyy 2 A(y)+x / dy 

J0 Jx 



A(y)F(y) 



(24) 



From this last equation it is easy to see (e.g., by means of the change of variables y = zx) that if A 
and F are IR finite, then L(0) = 0, as mentioned before. 

At this point one can substitute into the equations above the available lattice data on the gluon 
and ghost propagators thus determining in an indirect way the functions u and L. The results 



2 It is easy to recognize, for example, by substituting into ( p2p tree-level expressions, that (q 2 ) and u(q 2 ) have 
the same leading dependence on the ultraviolet cutoff A, namely F^ l (q 2 ) = u vv (q 2 ) = 9g 2 /64K 2 \og(A 2 /q 2 ), while 
L(q 2 ) is finite (independent of A). 

3 The inclusion of the g^ v term is absolutely essential for the self-consistency of the entire renormalization pro- 
cedure. To be sure, the g^ v term appears naturally, given, for example, the form of the gluon propagator background 
quantum identity m 7B. 
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q 2 [GeV 2 ] 



Figure 5: The renormalization-group invariant product d(q 2 ) obtained combining the lattice results for the 
gluon propagator and our solutions for the function u(q 2 ) = G(q 2 ) according to Eq. (25) 



are shown in Fig. from which we explicitly see the ji -dependence of the KO function and in 
particular of the KO parameter, as well as the vanishing in the deep IR of the L(q 2 ) form factor. 
Notice that the KO function saturates in the deep IR at the value 4 u(0) ~ —0.6. 

Now, it is well-known that the product g 2 A(q 2 ) is a RG-invariant quantity, representing the 



non-Abelian generalization of the QED quantity e 2 A(q 2 ). Then using Eq. (15) we find that the 
product 

d[q) [\+u(q 2 ^ 2 )} 2 ' ^ 

is RG-invariant. But this then shows definitively that (z) the KO function cannot be an RG-invariant 
combination since (z'z) its /I dependence must be such that it cancels the pL -dependence of the 
numerator, as it is explicitly shown in Fig. g. 

The jx -dependence of the KO parameter is displayed in the left panel of Fig. ||; the observed }X- 
dependence is really sizable when contrasted with the impressive absence of any such dependence 



displayed by the genuinely RG-invariant quantity given in Eq. (25) which was computed using 
exactly the same sets of lattice data (Fig. ||). On the right panel of the same figure we show finally 
a comparison between our indirect determination of the KO function and the direct one obtained 
in [Q] where the function @ was studied by means of Monte Carlo averages: evidently the two 
curves compare rather well. 



Let us finally take a closer look at the background-quantum identity (|15|); due to the central 
equality of Eq. (^H), we have that in the Landau gauge 



4 A value for the KO parameter of —2/3 has been in fact predicted in | p3| | by studying how the presence of the 
Gribov horizon affects the KO criterion. However since there is a residual dependence on /i 2 of the KO parameter we 
consider this a coincidence due to the choice of the renormalization point rather than a fundamental prediction of the 
theory. 
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6 8 10 12 14 16 18 20 0,01 0,1 1 10 

n 2 [GeV z ] q 2 [GeV 2 ] 

Figure 6: Left panel: The dependence of the KO parameter u on the renormalization point ji; the red solid 
line corresponds to a fit of a phase transition type — u = a{p} — b) c , with parameters a = 0.633, b = 3.57, 
c = 0.025. Right panel: The KO function, —u(q 2 ), obtained from the solution of Eq. (24) (solid black line) 
compared to the lattice data of at jU = 4GeV 

Interestingly enough, this simple formula expresses the KO function in terms of two gluon prop- 
agators calculated in the Landau gauge of two very distinct gauge-fixing schemes, with no direct 
reference to the ghost sector of the theory. This observation opens up the possibility of deduc- 
ing the structure of the KO function using an entirely different, and completely novel, approach. 
Specifically, one may envisage a lattice simulation 5 of A; then, u(q 2 ) may be obtained from ( |26| ) by 
simply forming the ratio of the two gluon propagators. Given that A(0) is found to be finite on the 
lattice [JT], |2[] , it is clear that, in order for the standard KO criterion to be satisfied (i.e., u(0) = — 1), 
A must diverge in the IR. Needless to say, we consider such a scenario highly unlikely. What is 
far more likely to happen, in our opinion, is to find a perfectly finite and well-behaved A, which in 



«GO = 0.388 and ;i = 3.0 GeV 
a(n z ) = 0.330 and p. = 3.6 GeV 
q(n') ; 0.295 and g ; 4.3 GeV 




q[GeVl 



Figure 7: The gluon propagator A(q 2 ) of the BFM, renormalized at three different points: ji = 3.0GeV 
(black curve), fl — 3.6GeV (red curve) and ji = 4.3GeV (green curve). 



5 The lattice formulation of the background field method has been presented in |g5p; interestingly enough, it was 
carried out in the Feynman gauge, which is the privileged gauge from the point of view of the pinch technique um. 
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the deep IR will be about an order of magnitude larger than A(0), furnishing a value «(0) ~ —0.6, 
namely what we have found in our analysis. In fact, one may turn the argument around: combining 
the results of this article with the lattice data for A [Jxj, one may use (26) to predict the outcome 



of the lattice simulation for A; our prediction for the case of SU (3) is shown in Fig. [7|. 
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